Yet another surprise in the problem of classical diamagnetism 
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The well known Bohr-van Leeuwen Theorem states that the orbital diamagnetism of classical 
charged particles is identically zero in equilibrium. However, results based on real space-time ap- 
proach using the classical Langevin equation predicts non-zero diamagnetism for classical unbounded 
(finite or infinite) systems. Here we show that the recently discovered Fluctuation Theorems, namely, 
the Jarzynski Equality or the Crooks Fluctuation Theorem surprisingly predict a free energy that 
depends on magnetic field as well as on the friction coefficient, in outright contradiction to the 
canonical equilibrium results. However, in the cases where the Langevin approach is consistent with 
the equilibrium results, the Fluctuation Theorems lead to results in conformity with equilibrium 
statistical mechanics. The latter is demonstrated analytically through a simple example that has 
been discussed recently. 
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I. INTRODUCTION 

The well known Bohr-van Leeuwen Theorem (BvL) 
GL H II II H IS S H forbids the presence of orbital dia- 
magnetism in classical equilibrium systems. The essen- 
tial point of the proof of this theorem is that the mag- 
netic field B enters the particle Hamiltonian through the 
replacement of the particle momenta p by p + e ^ , 
where A(r) is the associated vector potential and — e is 
the charge of the particle. Since the partition function 
involves integration of the momenta over the entire mo- 
mentum space, the origin of p can be trivially shifted by 
eA J r ) , and as a result, A(r) disappears from the partition 
function. This in turn implies that the free energy un- 
dergoes no change in the presence of a magnetic field and 
hence gives zero orbital magnetic moment. This result is 
rather surprising, given the fact that each particle must 
trace a cyclotron orbit in the presence of a magnetic field 
and, therefore, contribute to the diamagnetic moment. 
This was resolved by noting that the skipping orbits of 
the electron at the boundary generate paramagnetic mo- 
ment equal and opposite to that due to a carrier in the 
bulk 0,1, 1 IH. Thus the bulk diamagnetic contribu- 
tion is exactly cancelled by the boundary (surface) con- 
tribution, leading to total absence of orbital magnetism 
in classical equilibrium systems. The canonical statis- 
tical mechanical treatment, however, makes no explicit 
reference to such boundary effects. 

Now, let us consider a finite/infinite classical system 
where the particle does not hit a geometrical boundary 
all along its motion. In such a situation, classical diamag- 
netism is expected as the skipping trajectories carrying 
paramagnetic current along the boundary are absent [j| . 



This subtle role of the boundary has been revisited by 
Kumar and Kumar [9| by considering the motion of a 
charged particle which is constrained to move on the sur- 
face of a sphere, i.e., on a finite but unbounded system. 
The surface of a sphere has no boundary, and to the pleas- 
ant surprise of the authors, they did find non-zero classi- 
cal orbital diamagnetic moment by following the space- 
time approach. This effect has been attributed to the 
dynamical correlation induced by Lorentz force between 
velocity and transverse acceleration when the problem is 
treated as per the Einsteinian approach, i.e., in this case, 
via the Langevin dynamics [Tol . |Tl| . Such subtle dynami- 
cal correlations are presumably not captured by the clas- 
sical Gibbsian statistical mechanics based on equilibrium 
partition function. 

In our present work, we explore this system further us- 
ing the recently discovered fluctuation Theorems (FTs), 
namely, the Jarzynski Equality (JE) and the Crooks' 
Fluctuation Theorem (CFT) These FTs address 

the calculation of equilibrium free energy difference AF 
between two thermodynamic states derivable from irre- 
versible (nonequilibrium) trajectories. We come across 
other intriguing consequences. If the system is driven out 
of equilibrium by perturbing its Hamiltonian (H\) by an 
externally controlled time-dependent protocol A(i), the 
thermodynamic work done on the system is given by [l2j] 



W 



■dH , 



(1) 



over a phase space trajectory, where r is the time through 
which the system is driven. A(0) = A and A(r) = B are 
the thermodynamic parameters of the system. The JE 
states 
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(e-P w ) = e -P AF , 



(2) 



where AF = Fb ~ Fa is the free energy difference be- 
tween the equilibrium states corresponding to the ther- 
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modynamic parameters B and A, and the angular brack- 
ets denote average taken over different realizations for 
fixed protocol X(t). In eq.©, /3 = l/k B T, T being the 
temperature of the medium and ks is the Boltzmann 
constant. Initially the system is in equilibrium state de- 
termined by the parameter A(0) = A. The work done 
W during each repetition of the protocol is a random 
variable which depends on the initial microstate and on 
the microscopic trajectory followed by the system. The 
JE acts as a bridge between the statistical mechanics of 
equilibrium and noncquilibrium systems and has been 
used experimentally [lj] to calculate free energy differ- 
ences between thermodynamic states. The CFT predicts 
a symmetry relation between work fluctuations associ- 
ated with the forward and the reverse processes under- 
gone by the system. This theorem asserts that 



which in polar coordinates reduces to 



PfiW) 

Pr(-W) 



o 0(W-AF) 



(3) 



where PfiW) and P r (W) denote distributions of work 
values for the forward and its time-reversed process. Dur- 
ing the forward process, initially the system is in equi- 
librium with parameter A. During the reverse process, 
the system is initially in equilibrium with parameter B 
and the protocol is changed from Ab to Xa over a time 
r in a time reversed manner (X(t) — A(r — t)) and in our 
present problem, magnetic field also has to be reversed in 
sign Q. From equation ©, it is clear that the two dis- 
tributions cross at W = AF, thus giving a prescription 
to calculate AF. 

In the present work, we show that the case of a charged 
particle moving on a sphere leads to free energy of the sys- 
tem which depends on the magnetic field and on the dis- 
sipative coefficient, which is inconsistent with the predic- 
tion of canonical equilibrium statistical mechanics. The 
same system gives orbital diamagnetism when calculated 
via the space-time approach, again in contradiction with 
the equilibrium statistical mechanics p]. For the re- 
cently studied case of a particle moving on a ring [l5l |. 
the Langevin approach predicts zero orbital magnetism, 
just as in the present treatment. Thus, in this case, the 
free energy obtained by using FTs is consistent with the 
canonical equilibrium statistical mechanics. 



II. CHARGED PARTICLE ON THE SURFACE 
OF A SPHERE 

We take up the model proposed in [9], which consists 
of a Brownian particle of charge — e constrained to move 
on the surface of a sphere of radius a, but now with a 
time-dependent magnetic field B(i) in the z direction. 

The Hamiltonian of the system in the absence of heat 
bath is given by: 



H = 



1 

2m 



eA(r,t) 



(4) 



H 



1 

2m 



po_ 

a 



eA e (t) 



zA<f,{t) 



a sin 9 



(5) 

In a symmetric gauge, Ag — and A^ = 
(l/2)aB(t) sm9. In presence of the heat bath, the dy- 
namics of the particle is described by the Langevin equa- 
tion % 



dv e, „, „ e dBit) 
m— = — vxBi -IV rx — ±± 

dt c y " 2c V dt 

+V27T£(i), 



(6) 



where m is the particle mass and T is the friction coef- 
ficient. £(t) is a Gaussian white noise with the prop- 
erties = and (&(t)&(f )) = S M S(t - t'). The 

first term on the right hand side is the Lorentz force. 
If the magnetic field varies with time, it also produces 
an electric field E, hence the presence of the force term 
-eE = -(e/2c)(rx ^§(t)) ineq©. This is an additional 
element of physics not present in reference Q . 

Switching over to the spherical polar coordinates 0], 
eq. ^ assumes the following form in terms of dimen- 
sionless variables: 



sin f cos ( 



> sin V cos ( 



07 



(7a) 



) sm f 



cos ( 



au> c (B(t)) ■ „ ab ■, . . „ 

— v y!1 e cos 6 H Bit) sin 9 

c c 



07 



> sin 9 + ^Jrj , 



(7b) 



In the above equations, the dots represent differentia- 
tion with respect to the dimensionless time r = (c/a)t. 
Here 7 = T/m, ui c (B(t)) = eB(t)/mc, b = e/(2mc) and 
r\ = 2Taj/mc 3 . 



III. RESULTS AND DISCUSSIONS 

First we consider the case of static magnetic field B of 
magnitude B in the z direction. The ensemble averaged 
orbital magnetic moment which by symmetry is also in 
the z direction is given by 



(M(i)) 



ea 



(8) 



where (• • •) denote ensemble average over different real- 
izations of the stochastic process. 

Following the same procedure as in 0, we have cal- 
culated the equilibrium magnetic moment by double av- 
eraging first over a large observation time and then over 
the ensemble: 



3 



M eg = ((M(t))) = - f dt(M(t)) (9) 

T JO 

as t — > oo. For a numerical check, we have obtained the 
same results as figures 2 and 3 of @. Throughout our 
analysis, we have used dimensionless variables, e, c, m 
and a are all taken to be unity. 
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FIG. 1: Plots of magnetic moment M eq versus magnetic field 
B for different 7 and for a given temperature T — 1 . The 
different plots are for 7=1, 1.5 and 2, as mentioned in the 
figure. 
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FIG. 2: Plots of M eq versus B for different T and for a given 
friction coefficient 7=1. We have taken different plots for 
T=l, 1.5 and 2. 

The thermodynamic work done by the external time- 
dependent magnetic field on the system up to time t is 
given by 

W{t) = f ^§dt' = ^- fdt' j>(t')sm 2 6(t')B(t'). 

(10) 

In our case, B(t) acts as the external protocol \{t). 
The Langevin equations are solved numerically by using 
the Euler method of integration with time step At = 




FIG. 3: Plots of Meq as a function of T for 7 = 1 and for 
three different values of the external magnetic field: B=3, 5 
and 7. 
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FIG. 4: Plots of M eq as a function of the friction coefficient 
7, for 4 different values of B: B=3, 5, 7 and 10, with T=l. 
Note that the axis for 7 starts from 0.8. In the inset we have 
plotted the curves M eq versus 7 for B=12 and 15. 



0.01. Same boundary conditions and numerical proce- 
dure is carried out as in Q. 

In figure[lj we have plotted the dimensionless magnetic 
moment M eq {= -^-) versus the magnetic field in dimen- 
sionless units B{= ^f) for different values of the friction 
coefficient 7(= — ), as mentioned in the figure. At each 
point, the signature of M eq is opposite to that of B, pro- 
viding clear evidence of diamagnetism. Initially, M eq in- 
creases with B (linear response) and after showing a peak 
at high fields, it approaches zero. At high fields, it is ex- 
pected that the radius of the cyclotron orbits will tend 
towards zero, and hence naturally the magnetic moment 
also vanishes. With increase in the friction coefficient 7, 
the peak shifts towards higher magnitudes of magnetic 
field. It should be noted that this behaviour is qualita- 
tively consistent with the exact result obtained for the 
orbital magnetic moment M.1& for a charged particle in 
a two-dimensional plane in the absence of a boundary, 
following the real space-time approach (see eq.(8) of [f|). 
The expression for orbital magnetic moment M.%s. is given 
by 
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FIG. 5: Plots of AF versus the final value of the magnetic 
field B(t) for 7 = 1 and 7 = 2. The protocol used is a 
ramp, B = Bot/r, for a time of observation r = 2000, with 
the temperature fixed at T=l. The inset shows the variation 
of AF as a function of the friction coefficient 7, with the 
parameters T=l, B(t) — £?o=10. 










„ P(-W) 



T=1,T=2000, B(T)=15 




FIG. 6: Determination of AF using the CFT. (a) Plots of 
Pf(W) and P r (-W) at B(r) = B = 10, which cross at W(= 
AF) = 0.22. (b) Plots of P f {W) and P r (-W) at B{t) = 
B = 15, which cross at W{= AF) = 0.65. 



M M = - — 



2c V7 



Tuj c 



(11) 



Here, tu c — eB/mc, where B is the magnitude of the 
static magnetic field. Compared to the analysis in [§], 
we have gone beyond the linear response regime. 

In figure [5] we have plotted the magnitude of M eq as a 
function of B for different values of temperature T. From 
figures Q] and [2j it can be inferred that the magnetic mo- 
ment can be monotonic or non-monotonic in T and 7, 
depending on the whether the values of B lie within the 
linear response regime or beyond. To this end, in figures[3] 
and[H we have plotted the equilibrium magnetic moment 
as a function of temperature T and friction coefficient 7 
respectively, for various values of B. The magnetic mo- 
ment is zero at T = as well as at T = 00. It exhibits a 
minimum in the intermediate range of temperature. This 
minimum shifts towards lower temperature with the in- 
crease in B. It should be noted that for larger temper- 
atures, a higher number of realizations are required to 
generate more accurate data points. 

In figured! we notice that in the parameter range that 
we have considered, the equilibrium magnetic moment 
decreases monotonically with friction coefficient. For 
large 7, the particle motion gets impeded by the medium 
and as expected, M eq — > as 7—5-00. As 7 — > 0, there is 
a saturation in the value of magnetic moment, which de- 
pends on the value of the parameters B and T. This we 
have not shown in the figure. It is evident from figure [1] 
that for large values of B (B > 10), dependence of M eq on 
7 is non-monotonic. This is shown in the inset where M eq 
is plotted as a function of 7 for B — 12 and for B = 15. 
It is observed that the dip in M eq shifts towards higher 
7 for higher value of B. For small friction coefficients, 
the saturation value is very small (for large B) and it 
requires a much larger number of realizations to achieve 
reliable results. The details of these results will be pub- 
lished elsewhere. Our results clearly indicate that the 
temperature and the friction dependence of the classical 
magnetic moment obtained via real space-time approach 
are qualitatively different for an infinite unbounded sys- 
tem (eq. (HU) from that for a finite unbounded system 
considered here. From eq. (|11[) we can readily infer that 
the dependence of on temperature T and on friction 
coefficient 7 is monotonic. 

Having shown that the space-time approach leads to 
a finite diamagnetic moment in contrast to its absence 
in canonical equilibrium, we can now turn to the calcu- 
lation of free energy differences for the same problem 
using the FTs. We subject the system to the time- 
dependent magnetic field (protocol) in the form of a 
ramp, B(t) — Bot/r, where r denotes the total time 
of observation. We use the ramp with an observation 
time r = 2000. The final value of magnetic field is 
B(t) — Bq. To calculate the free energy difference, 
AF = F(B ) - F(0), we have used the JE (eq ©). To 
calculate AF numerically , we have generated 10 4 real- 
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izations of the process, making sure that the system is 
initially in canonical equilibrium in the absence of mag- 
netic field (5(0) = 0). The results for AF are plotted 
as a function of B(t) in figure [5] for two values of 7. 
All physical parameters are in dimensionless units and 
are as mentioned in the figure. Surprisingly, we notice 
that AF depends on the magnetic field B(t). This is 
in sharp contrast to the equilibrium result, namely, AF 
should be identically zero. To our knowledge, this is the 
first example wherein the Fluctuation Theorem fails to 
reproduce the result obtained from equilibrium statisti- 
cal mechanics. This is yet another surprise in the field of 
classical diamagnetism. Moreover, AF depends on the 
type of protocol. The dependence of AF on the friction 
coefficient is shown in the inset of figure \5\ In classi- 
cal equilibrium, it should be noted that the free energy 
does not depend on friction coefficient. From this free 
energy, one can get moment by calculating the derivative 
of the obtained free energy with respect to B. However, 
the magnetic moment thus obtained does not agree with 
that obtained through the simulation of the Langevin 
equations. This we have verified separately. 

In figure [5] (a) and (b), we have plotted Pf(W) and 
P r (—W) as a function of W for the same protocol ending 
with two different values of the magnetic field (B(t) = 10 
and 15). The crossing point of Pf(W) and P r (— W), ac- 
cording to the CFT, gives the value of AF, which we 
have found to be equal to 0.22 for Bo — 10 and 0.65 
for Bq — 15, which are in turn equal to the obtained 
values using the JE, namely, 0.22 and 0.65 respectively, 
within our numerical accuracy. Thus we have shown that 
a charged particle on a sphere exhibits finite diamagnetic 
moment and magnetic field dependent free energy calcu- 
lated via real space-time approach and the Fluctuation 
Theorems respectively. As mentioned earlier, these re- 
sults contradict equilibrium statistical mechanics. 



IV. CHARGED PARTICLE ON A RING 

Now we turn to a simpler problem of a charged parti- 
cle moving on a ring in a magnetic field perpendicular to 
the plane of the ring, i.e., in the z direction. This prob- 
lem has been studied recently [l5| in connection with 
the BvL for a particle motion in a finite but unbounded 
space, where it was shown that this system analyzed via 
the Langevin dynamics does not exhibit orbital diamag- 
netism, consistent with equilibrium statistical mechanics. 
It is not surprising as the equation of motion for the rel- 
evant dynamical variable, namely the azimuthal angle 
<f>, does not depend on the strength of the static mag- 
netic field. Hence, the magnetic field has no effect on the 
motion of a particle constrained to move in a circle of 
fixed radius a. We analyze the same problem, however in 
the presence of time-dependent magnetic field (protocol) , 
within the framework of Jarzynski Equality to obtain the 
free energy dependence on magnetic field in this case. To 
this end, the Hamiltonian of the system is given by 



H 



J_ fp± , eA^t) 
2m \ a c 



(12) 



where, for a magnetic field in the z direction, A^it) — 
(a/2)B(t). The corresponding Langevin equation for the 
relevant variable <p is given by 



TOa = _r a ^+— B(t)+V2IT (13) 
The above equation can be written in a compact form 



(14) 



with 7 = — , A = tt~ — and 77 = ^t-. In this section, the 

' m ' 2mc ' raa A ' 

dots represent differentiation with respect to real time t. 
It may be noted that if the magnetic field is independent 
of time, i.e., B = 0, then the field has no effect on the 4> 
variable, as can be seen from equation (|14[). The thermo- 
dynamic work W, using equation (fT]) and (TT2"j) . is given 
by 



,2 ft 



W(t) = I ^dt' = — <t>(t')B{t') dt'. (15) 



r * dH 
lo ~W 

The formal solution for <fi is given by 



ea 

~2c J Q 



4>{t) = <K0)e" 



dt' [XB(t') + 07 £*(*')]• 
(16) 



Substituting this solution in eq (fTSf for W, we get 



W{t) = g / dt'B(t')[(j>(0)e 



-if 



{\B(t") 



+V^0(O}e 7t " dt"], 



(17) 



where g = ea 2 /2c. 

Since the expression for W in the above equation is 
linear in the Gaussian stochastic variable £,^(t), W it- 
self follows a Gaussian distribution. To obtain P(W), 
we simply need to evaluate the average work (W) and 
the variance = (W 2 ) — (W) 2 . the full probability 
distribution P{W) is given by 



P(W) 



y/27 



: exp 



(W — (W)) 2 



(18) 



Averaging eq. (jTTJ) over random realizations of £<f,(t), 
and noting that (^(t)) = 0, we get for average work done 
till time r: 



(W)=gX dt'Bit'^-^' dt" '. (19) 

Jo Jo 
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Again using eq. (fT7|) and (TT9]) . after tedious but 
raightforward alg 
obtained and is given by 



straightforward algebra, the variance can be readily 



independent of magnetic field. We have obtained this 
exact result which is independent of the functional form 
the protocol B(t). 



°w = 7TV I dt' B{t') [ dt 1 B(t 1 )e-Tl*'-* I l. (20) 
2 1 Jo Jo 

In arriving at the above expression, we have used the 
fact that the variance of the initial equilibrium distribu- 
tion of angular velocity </>(0) is given by (</> 2 (0)) = ^7 = 
\g 2 T). Comparison between (fTO)) and (|2"0|) gives the result 

a 2 w = 2T{W), (21) 

a fluctuation-dissipation relation. 
Using eq. (fTSJ) and (|2"T]). we get 

(e~ fiW ) = 1, (22) 

which, according to the JE, implies AF = F(B(t)) — 
F(B(0)) = 0, where B{0) and B(t) are the values of the 
magnetic field at the initial and final times of the proto- 
col respectively. The magnitudes of B(0) and B(t) = B 
can take any value. Thus, AF = implies that the free 
energy is independent of the magnetic field, the result 
being consistent with equilibrium statistical mechanics. 
It is interesting to note that the averaged work (W) (eq. 
(fni)0 and its variance cr 2 ^ (eq. (f2"0"]0 depend on the func- 
tional form of B(t) and on the magnetic fields at the end 
points of the observation time and yet (exp(— (3W)) is 



V. CONCLUSION 

In conclusion, whenever the real space-time approach 
for a charged particle in the presence of a magnetic 
field predicts a finite diamagnetic moment, the Fluctu- 
ation Theorems too fail to reproduce results consistent 
with equilibrium statistical mechanics. These conclu- 
sions have also been supported by the results for the 
motion of a charged particle in a two-dimensional plane 
in the absence of boundary 0, [l6| . In cases where real 
space-time approach to diamagnetism is not in conflict 
with the equilibrium statistical mechanics, an example 
being a charged particle on a ring or in the presence of 
a confining boundary Q , the Fluctuation Theorems lead 
to results consistent with equilibrium statistical mechan- 
ics. Only experiments can resolve whether really orbital 
diamagnetism exists in classical equilibrium systems (like 
charged particle on the surface of a sphere). 
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